Fractional Chern Insulators in Topological Flat bands with Higher Chern Number 
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Lattice models forming bands with higher Chern number offer an intriguing possibility for new phases of 
matter with no analogue in continuum Landau levels. Here, we establish the existence of a number of new bulk 
insulating states at fractional filling in flat bands with Chern number C — N > 1 , forming in a recently proposed 
pyrochlore model with strong spin-orbit coupling. In particular, we find compelling evidence for a series of 
stable states at v — 1/(2N + 1) for fermions as well as bosonic states at v = 1/(N + 1). By examining 
the topological ground state degeneracies and the excitation structure as well as the entanglement spectrum, 
we conclude that these states are Abelian. We also explicitly demonstrate that these states are nevertheless 
qualitatively different from conventional quantum Hall (multilayer) states due to the novel properties of the 
underlying band structure. 



PACS numbers: 73.43.Cd, 71.10.Fd, 73.21.Ac 

Introduction. — The recent discovery of nearly fiat bands 
with unit Chern number C = 1 in itinerant lattice systems [1- 
3] has spurred plenty of theoretical excitement [4-20] as these 
bands may harbor lattice analogues of fractional quantum Hall 
(FQH) states that do not require an external magnetic field and 
may potentially persist at very high temperatures. 

While flat C = 1 bands are very similar to continuum Lan- 
dau levels, lattice models can harbor bands with higher Chern 
number and may therefore host qualitatively new phases of 
matter with no analogue in the continuum. Two independent 
groups have very recently shown that fiat bands with arbitrary 
higher Chern number can in fact be systematically created in 
multilayer systems [21, 22] assuming only short-range hop- 
ping (see also Ref. [23] for a nice earlier construction of a flat 
C = 2 band) and, as such, provide promising platforms for 
new phenomena [21-25]. 

In this work we study the crucial effect of interactions in the 
original pyrochlore (multi-layer kagome) lattice model pro- 
posed for flat bands with arbitrary Chern number C — N [21]. 
While a very recent study [24] focusing on interactions in a 
fairly flat C = 2 band on the triangular lattice reported an in- 
triguing v = 1/3 bosonic state, we find compelling evidence 
for a whole series of incompressible fractional Chern insula- 
tor (FCI) states, both for fermions at v = 1/(2N + 1) and for 
bosons at v = 1/(N + 1), N = 2, 3, 4. We also demonstrate 
that the states which we discover are, despite a number of sim- 
ilarities, qualitatively different from conventional (single and 
multilayer) FQH states. 

Setup. — We focus on the model describing Rashba spin- 
orbit coupled particles on pyrochlore slabs including N 
kagome layers (from /Ci to /Cjy) introduced in Ref. [21]. The 
single-particle Hamiltonian in real space is 

H = ]T tijcl^ + i ^2 A '.' :K ',' X R u) ' OafiCiaCjP- 
i,j,(r i,j,a,(3 

Although next-nearest hopping is needed to get very flat C = 



N bands, we consider only nearest neighbor hopping tij — 
t\ = — 1 and Xij — Ai (Ai = 1.1 for fermions and Ai = 
0.9 for bosons) within each kagome layer, as well as tij — 
t± [26] when involving the triangular layers for simplicity. 
In momentum space, there are AN — 1 bands assuming spin 
polarization, and there is a flat band with Chern number C — 
N for suitable choices of hopping parameters [21]. 

As is customary we take the flat band limit and consider the 
case when the flat band is partially filled by interacting parti- 
cles with an interaction Hamiltonian H mt . We diagonalize if; nt 
projected to the flat band for a finite system with N\ x N 2 unit 
cells. 

Fermion energetics. — We consider fermions with a near- 
est neighbor repulsion H m = Y2u j) n i n j [27] and begin by 
considering a bilayer kagome system at filling v = 1/5 in 
the C = 2 band. For each system size that we study, there 
are five quasidegenerate ground states at the bottom of the en- 
ergy spectrum. These lowest states are separated from other 
excited states by a clear energy gap, which is a necessary 
condition for the v = 1/5 fermionic FCI state [Fig. 1(a)]. 
The energy gap is always significantly larger than the ground 
state splitting for various system sizes, and a finite-size scal- 
ing analysis of the energy gap shows that it is very likely to 
survive in the thermodynamic limit [Fig. 1(b)] [28]. To cor- 
roborate that the ground states are topologically nontrivial, we 
also check the spectral flow under twisted boundary condi- 
tions, which amounts to inserting magnetic flux through the 
cycles of the system (i.e., through the handles of the torus). 
For a many-body state vp, the twisted boundary condition in 
the x{y) direction is ^(r^ + -/V 1 ( 2 )a 1 ( 2 )) = exp(i$)^(rj), 
where $ is the boundary phase and a 1( - 2 ) is the lattice vector. 
By calculating the spectral flow for some system sizes where 
the ground states are in different (Ki^K?) sectors, we find 
that when $ changes from to 5 x 2-7T, the five ground states 
evolve into each other, being always separated from excited 
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Figure 1. (Color online) Results for the v = 1/5, C = 2 fermionic 
FCI state in a bilayer kagome system with Ni — 5 and N2 = N e . (a) 
The low-lying energy spectrum for N e = 4, N e — 6, and N e = 8. 
(b) The finite-size scaling analysis for both energy gap and ground 
state splitting, (c) The x-direction spectral flow for N e = 7. (d) The 
quasihole excitations for N e = 4, Ni = 3, and N2 = 7 (one hole 
is added, 21 states below the gap), (e) The quasihole excitations for 
N e = 5, Ni = 3, and N2 = 9 (two holes are added, 81 states below 
the gap), (f) The quasihole excitations for N e = 5, Ni = 5, and 
N2 — 6 (five holes are added, 756 states below the gap). 



Figure 2. (Color online) Results for the v = 1/7, C — 3 fermionic 
FCI state in a three-layer kagome system with N 1 — N e and N2 = 
7. (a) The low-lying energy spectrum for N e = 5, N e = 6, and 
N e — 7. (b) The finite-size scaling analysis for both energy gap and 
ground state splitting, (c) The y-direction spectral flow for 7V e = 5. 
(d) The quasihole excitations for N e = 5, iVi = 6, and N2 = 6 (one 
hole is added, 36 states below the gap), (e) The quasihole excitations 
for N e — 4, Ni = 5, and N2 = 6 (two holes are added, 75 states 
below the gap), (f) The quasihole excitations for N e = 5, Ni = 6, 
and N2 — 7 (seven holes are added, 2772 states below the gap). 



states by a gap, and finally return to the initial configuration 
[Fig. 1(c)]. The behavior of the spectral flow indicates that 
the Hall conductance is an = f %- [29], which we have also 
confirmed by calculating the many-body Chern number. 

The quasihole excitations, which carry fractional charge 
[30] and statistics [31-33], are one of the most important char- 
acteristics for FQH states in Landau levels. Here, we first in- 
vestigate the quasihole excitations of our C = 2, v = 1/5 
fermionic FCI state. By keeping N e fixed and changing iVi 
and/or N2, we can add one hole and more holes into the sys- 
tem. A clear gap that separates the low-lying states from high- 
excited states exists in the quasihole excitation spectrum. The 
total number of states below the gap is identical to that for 
the v = 1/5 fermionic Laughlin state in the Landau level on 
the torus. Moreover, the ground state momenta and counting 
of quasiholes can be predicted by an exclusion rule known 
from the thin-torus limit of the FQH system [34]: By fold- 
ing the two-dimensional momenta (ki, k^) to one dimension 
as fciD = ki + Nxk2 [5], this rule implies that there are no 
more than p particles in q consecutive orbitals at, and slightly 
below, v = p/q. This shows that the quantum dimension of 
the quasiholes is d q h = 1, which is a hallmark of Abelian 
statistics [35]. 

In the three-layer kagome system, we find similar evidence 
for a v = 1/7, C — 3 fermionic FCI state (see Fig. 2 for a 
summary of these results). The gap scaling is again convinc- 
ing albeit with a gap being roughly 20% of the gap observed at 



v = 1/5 in the C = 2 band. A clear gap in the quasihole exci- 
tation spectrum of the v = 1/7 state can be identified, and the 
total number of states below the gap is identical to that for the 
v = 1/7 fermionic Laughlin state in the Landau level on the 
torus. In the four-layer kagome system, we find some clues of 
v = 1/9, C = 4 fermionic FCI state, again clearly resolved, 
albeit with a yet smaller gap [36]. 

Boson energetics. — We now turn our attention to bosons 
with on-site repulsion i7; nt = X)i n i( n « — 1) m tne bilayer, 
three-layer, and four-layer kagome systems. In the bilayer 
kagome system, we find af = 1/3 bosonic FCI state [36], 
presumably corresponding to the same phase as the v = 1/3 
recently observed in a C — 2 band on the triangular lattice 
[24]. In the three-layer kagome system, we find convincing 
evidence for a new v = 1/4 bosonic FCI state in the C = 3 
band as displayed in Fig. 3 [37]. We also find an interesting 
v = 1/5, C = 4 bosonic FCI state in the four-layer kagome 
system [36]. Note that, in contrast to C = 1 systems, bosonic 
FCI states can form at a Laughlin fraction with an odd de- 
nominator when C > 1. In each case, a clear gap exists 
above the ground state manifold, as well as in the quasihole 
excitation spectrum. Again, the counting of quasiholes can be 
predicted by the thin-torus-like exclusion rule. While the rule 
for predicting ground state momenta of the boson states must 
be slightly modified [38], this is analogous with the case in 
C = 1 bands. 

Entanglement spectra. — To further corroborate our find- 
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Figure 3. (Color online) Evidence for the v = 1/4, (7 = 3 bosonic 
FCI state in a three-layer kagome system with N\ = 4 and N^ = 
Nb- (a) The low-lying energy spectrum for Nb = 4, Nb = 6, and 
Nb = 8. (b) The finite-size scaling analysis for both energy gap and 
ground state splitting, (c) The x-direction spectral flow for Nb = 7. 
(d) The quasihole excitations for Nb = 5, Ni = 3, and N2 = 7 (one 
hole is added, 21 states below the gap), (e) The quasihole excitations 
for Nb = 4, Ni = 3, and N2 = 6 (two holes are added, 45 states 
below the gap), (f) The quasihole excitations for Nb = 5, Ni = 4, 
and A^2 = 6 (four holes are added, 336 states below the gap). 





Ml 



21 28 35 42 49 
K +NK 



Figure 4. (Color online) Particle entanglement spectra probing the 
Na = 3 quasihole excitations of the (a) N e — 7, v = 1/5, C = 2 
state on the 5x7 lattice (2695 states below the gap) and (b) iV e = 
7, v = 1/7, C = 3 state on the 7 x 7 lattice (7105 states below the 
gap)- 



ing of FCI phases, we have investigated the particle entan- 
glement spectra [39, 40]. This provides an independent test 
of the excitation structure of the system and can be used to 
discard competing possibilities, such as charge density waves 
[5]. The results from the largest systems we have studied are 
displayed in Fig. 4. Here, we find that, although the band- 
width of the manifold of low-energy levels is relatively large 
(which appears to be the generic situation for FCI as well as 
FQH states), there is a clear entanglement gap separating these 
levels from generic ones. Strikingly, the number of low-lying 
levels exactly matches the pertinent quasihole in the corre- 
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Figure 5. (Color online) (a) Energy gap (solid line) and ground 
state splitting (dotted line) as a function of a at v = 1/7 for var- 
ious system sizes: N e — 4,Ni = 4,N 2 = 7 (blue square); 
N e = 5,Ni = 5,N 2 = 7 (red cross); iV e = 6, Ni = 6, N 2 = 7 
(green triangle); and N e = 7, Ni — 7, iVa = 7 (yellow circle), (b) 
The average overlap \ £Li |<*<(o)|*<(a = 0))|, where |*i(a)> 
is the i th degenerate ground state. 



sponding Abelian FQH states. For smaller systems we some- 
times observe deviations when Na is too large [36], but the 
overall picture conclusively rules out a charge density wave 
explanation (cf., e.g., Ref. [20]). 

Consequences of mutual interlayer entanglement. — It is 
tempting to associate the higher Chern bands with conven- 
tional multilayer FQH systems as these systems are topologi- 
cally equivalent (both have N chiral edge states, etc.). How- 
ever, the dynamics of all layers in the higher Chern bands is 
invariably entangled (independent of the interlayer tunneling 
strength t± ^ 0) due to the novel band structure, while the 
layers in the FQH system can in principle be independent. 
To probe the consequences of this, we now consider the case 
of variable interaction strength in the different kagome lay- 
ers with fermions at v = 1/7 filling in a three-layer kagome 
system as an example (Fig. 5). Here we weaken the in- 



teraction in the midlayer K-2 as H mt (a) = 



<M>e/c 2 ' H 



ij with a G [0, 1] 



(i,5) 



We find that, even when the interaction in the midlayer is 
vanishing (a = 1), the ground state degeneracy and energy 
gap are still finite, and the ground states have a large over- 
lap with those for a = 0. We also calculate the spectral 
flow, the quasihole excitation spectrum, and the entanglement 
spectrum, and those results confirm that the ground states for 
a = 1 are qualitatively the same as those for a = [36]. 
This is strikingly different compared to the FQH multilayer 
case where the gap would close (at least in the weak tunneling 
regime). 

Discussion. — In this work, we have provided evidence for 
a number of new strongly correlated states emerging in nearly 
flat bands with higher Chern number. These states are likely 
to be Abelian, as they share ground state degeneracies and ex- 
citation structure with Laughlin-like states at the same filling 
fractions (whenever they exist). In addition to the fact that 
the FCI states found here can be realized at anomalous filling 
fractions (for bosons) compared to single-layer FQH states, 
we have shown that these new states are also qualitatively dis- 
tinct from conventional multilayer FQH states by weakening 
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the interaction in the midlayer. 

We note that, while multilayer FQH states generically have 
higher ground state degeneracy than the states found here, 
it is well known that, e.g., symmetrizing simple Abelian 
multilayer states can have dramatic consequences [41], most 
saliently leading to non- Abelian phases with a reduced ground 
state degeneracy. Symmetrizing can also lead to Abelian sin- 
gle band FQH states with a yet smaller ground state degen- 
eracy, reminiscent of the FCI states discovered here. In this 
context, we note that there are a number of intriguing similar- 
ities between the states reported here, and composite fermion 
states for bosons at v — and fermions at v = „ p , . , such 

p+l 2p+l' 

as sharing the same Hall conductance, ground state degener- 
acy, and filling fraction denominator. Moreover, the compos- 
ite fermion states have the form of (anti)symmetrized multi- 
layer states, and such states have recently been confirmed to 
exist also in C = 1 Chern bands [10, 11]. This is suggestive 
of a flux attachment picture [42] also for the C > 1 states 
presented here. 

It is very likely that the C > 1 bands can also harbor many 
other incompressible phases, and an exhaustive investigation 
of the phase diagram remains a challenge for future works us- 
ing more sophisticated techniques including considering tilted 
samples [11] to obtain a more precise finite-size scaling con- 
trol. Nevertheless, the states presented here appear to be the 
most stable ones as indicated by preliminary scans of more 
generic filling fractions (for relatively small system sizes). 

Another open issue is whether our findings are relevant to 
experiments. One class of candidate systems is cold atom se- 
tups with artificial gauge fields [43^4-5], while another that 
is especially well suited for the present model is provided by 
conveniently grown spin-orbit coupled solid state materials, 
such as the pyrochlore iridates [46-50]. 

Note added. — Recently, a pre-print with closely related re- 
sults for bosons, including further C > 1 FCI states in the 
model studied here, appeared [38]. 
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SUPPLEMENTARY MATERIAL FOR "FRACTIONAL CHERN INSULATORS IN TOPOLOGICAL FLAT BANDS WITH 

HIGHER CHERN NUMBER" 



In this supplementary material for "Fractional Chern Insulators in Topological Flat Bands with Higher Chem Number", we 
show the energetics of the v = 1/9 fermionic, v = 1/7 fermionic (with weakened interaction in the midlayer), v = 1/3 bosonic 
and v = 1/5 bosonic FCI states in the multilayer kagome lattice. For several of these systems we also present data on the 
particle entanglement spectrum (PES). We also give a systematic comparison between the PES counting of the FCI states and 
the counting of quasihole excitations derived from a simple exclusion rule (which gives the counting associated with Laughlin 
states at filling fractions where these exist). 



Evidence for the v = 1/9, C = 4 fermionic FCI states in the four-layer kagome system 



Here we show some data about the existence of the v — 1/9 fermionic FCI states in the four-layer kagome system, including 
the energy spectra and spectral flow (Fig. 6). One can see that there are 9 nearly degenerate ground states in the energy spectra, 
which are separated from other excited states by an energy gap (about 2 x 1CP 4 ~ 3 x 10~ 4 ). The spectral flow under twisted 
boundary condition also indicates that the ground states are topologically non-trivial and robust. 
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Figure 6. (Color online) Results for the v = 1/9, C — 4 fermionic FCI states in the four-layer kagome system (Ai = 1.1). (a) The low-lying 
energy spectrum for N e = 5(iVi = 5, iV 2 = 9) and N e = 6(A r i = 6, N2 = 9). (b) The y-direction spectral flow for N e = 5(JVi = 5, iV 2 = 
9). 



Evidence for the v = 1/5, C = 4 bosonic FCI states in the four-layer kagome system 



We also find an interesting v — 1/5 bosonic FCI state in the four-layer kagome system (Fig. 7). This is a new and somewhat 
unexpected bosonic state at a Laughlin fraction with odd denominator. A clear gap exists above the ground state manifold, as 
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well as in the quasihole excitation spectrum. The total number of states below the gap is consistent with that for the v = 1/5 
fermionic Laughlin state in Landau level. The energy gap decreases in the last data point in Fig. 7(b). But given the systematics 
of our findings we still believe that there is an FCI also in this case and that the drop is accidental for this particular system size. 

We also calculate the PES of the v = 1/5 bosonic FCI states (see Fig. 8, in which we also show the PES of the v = 1/4 
bosonic FCI states as well for completeness). A clear entanglement gap can be observed in the PES. A detailed discussion of 
PES can be found in the last section of this Supplementary Material. 
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Figure 7. (Color online) Evidence for the v — 1/5, C = 4 bosonic FCI state in the four-layer kagome system with N2 = 5 and Ni = Nb 
(Ai = 0.9). (a) The low-lying energy spectrum for Nb = 5, Nb — 6, and Nb = 7. (b) The finite-size scaling analysis for both energy gap and 
ground state splitting, (c) The {/-direction spectral flow for Nb = 6. (d) Nb = 4, Ni = 3, and N2 — 7 (one hole is added, 21 states below 
the gap), (e) Nb = 5, Ni = 3, and ./V2 = 9 (two holes are added, 81 states below the gap), (f) Nb = 5, Ni = 6, and N2 = 5 (five holes are 
added, 756 states below the gap). 



Evidence for the v = 1/3, C = 2 bosonic FCI states in the bilayer kagome system 



Bosonic v = 1/3 FCI states were very recently discovered in C = 2 band in a triangular lattice model [24]. We also observed 
such states in our C = 2 bilayer kagome system model. The ground state degeneracy and spectral flow are shown in Fig. 9. 



Spectral flow, quasihole excitations and entanglement spectrum for v = 1/7, C = 3 fermionic FCI states with weakened interaction 

To further confirm that the ground states are FCI states, we calculate the spectral flow, quasihole excitations and PES for 
v = 1/7 ground states with weakened interaction in the midlayer. We consider an extreme case in which the interaction in the 
midlayer is completely eliminated (a — 1). In Fig. 10, one can find that the spectral flow, quasihole excitation spectra and PES 
are qualitatively the same as those for a = 0. 
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Figure 8. (Color online) (a) The PES for the v = 1/5, C = 4 bosonic FCI state in the four-layer kagome system with Nt — 7, Ni — 7, and 
N 2 — 5 (2345 states below the gap), (b) The PES for the v = 1 /4, C — 3 bosonic FCI state in the three-layer kagome system with Nb = 7, 
iVi = 4, and N 2 = 7 (1428 states below the gap). Ai = 0.9. 
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Figure 9. (Color online) Results for the v — 1/3, C = 2 bosonic FCI states in the bilayer kagome system (Ai = 1) (a) The low-lying energy 
spectrum for N b = 4(iVi = 3, N 2 = 4), N b = 6(Ni = 3, N 2 = 6), and N b = 8(Ni = 4, N 2 = 6). (b) The x-direction spectral flow for 

^ = 4(^=3,7^2 = 4). 



Data of the PES counting for various filling factors and system sizes 



The particle entanglement spectrum (PES) is usually considered as a valuable tool to probe the excitation structure of the 
system and rule out competing states [5]. For a system possessing a d-fold (quasi)degenerate state we define the 

mixed density matrix of the system as p = \ Then we cut the total particles into two parts A and B with 7A^ 



(a) (b) (c) (d) 




Figure 10. (Color online) Results for the v = 1/7, C = 3 fermionic FCI states in the three-layer kagome system (Ai = 1.1) with weakened 
interaction H mt (a = 1), namely the interaction in the midlayer is eliminated completely, (a) The y-direction spectral flow for iV e = 5(7Vi = 
5, N 2 = 7). (b) The quasihole excitations for N e — 5, Ni — 6, and N 2 — 6 (one hole is added, 36 states below the gap), (c) The quasihole 
excitations for iV e = 4, iVi = 5, and ^2 = 6 (two holes are added, 75 states below the gap), (d) The PES for N e = 7, Ni = 7, and N 2 — 7 
(7105 states below the gap). 
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and Nb particles, respectively. The reduced density matrix of part A can be obtained by tracing out the particles belonging to 
part B, namely pa = Tr^p. We can label each eigenvalue of pa as e - ', where £ is just the PES level. As shown in some 
previous works, the PES of many FCI states in C = 1 band has a low-lying part separated from the higher levels by a clear 
entanglement gap and the number of levels in the low-lying part (PES counting) matches the quasihole excitations counting of 
the corresponding FQH states in Landau level. 

In this section, we show the PES counting for both fermionic and bosonic FCI states that we discover in higher C band and 
compare it with the quasihole excitations counting provided by the simple generalized exclusion rule, namely no more than one 
particle on m consecutive orbits (see Tab. I): 



kfm t N , _ ^HMMI m 
M ^ N ^~ „!(JV,-mrO! ' (1) 



where is the number of flux (orbits) and n is the number of particles. (This counting also holds for the v = 1/m Laughlin 
state on the torus.) 

After extracting PES from the FCI state at v = 1/m, we can get its counting Afp^ s (NA, Ni, N2). Then, we want to compare 
NpEs( N A,Ni,N 2 ) with Af^(N A ,N,p = NiN 2 ). In C = 1 band, they are found to always match [5]. Although the over 
all correspondence is striking in our higher C band simulations, we find some deviations. With the present data we cannot 
conclusively decide if these deviations have a physical meaning or if they are merely due to finite size effects. For fermionic 
v = 1/5 FCI states in C = 2 band, N^ S (N A , N U N 2 ) = N$(N A , = N 1 N 2 ). However, at other filling factors, in some 
cases where the system size is small and (or) Na is large, AT^^Na, N%, N 2 ) can be smaller than J\[^(Na, N$ = NiN 2 ). 
Despite of this derivation, our results clearly point to the Abelian nature of those FCI states and rule out the possibility of charge 
density wave states (for which the PES counting is expected to be significantly smaller [20]). 
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Table I. In this table, we give the PES counting of the FCI states for various filling factors and system sizes. For comparison, we also show 
the counting of the quasihole excitations provided by the simple generalized exclusion rule. 



fermions at v = 1/5 in C = 2 band 


svstem size ( l\f A/i y IVn') 


AP£*(Na Ni No) 


NT(Na Nj* — N No) 


iV e = 4, JV X x N 2 = 5 x 4 


JVa = 2 : 110 


N A = 2: 110 


iV e = 5, JVi x N 2 = 5 x 5 


N A = 2 : 200 


N A = 2 : 200 


N e = 6, Ni x N 2 = 5 x 6 


AT A = 2 : 315 


N A = 2 : 315 


A^a = 3 : 1360 


Na = 3 : 1360 


AT e = 7, JVi x N 2 = 5 x 7 


N A = 2 : 455 


N A = 2 : 455 


N A = 3 : 2695 


N A = 3 : 2695 


AT e = 8, JVi x N 2 = 5 x 8 


N A = 2 : 620 


N A = 2 : 620 


A^a = 3 : 4680 


N A = 3 : 4680 


N A = 4 : 17710 


AU = 4 : 17710 


fermions at i/ = 1/7 in C = 3 band 


system size (AT e , N 1 x A^) 


J^Ss(N A ,N u N a ) 


Af£(N A ,N <l> = N 1 N 2 ) 


jV e = 4, M x AT 2 = 4 x 7 


N A = 2: 210 


N A = 2: 210 


AT e = 5, iVi x AT 2 = 5 x 7 


A^a = 2 : 385 


A?a = 2 : 385 


AT e = 6, JVi x AT 2 = 6 x 7 


A^a = 2 : 609 


Na = 2 : 609 


A/a = 3 : 3248 


A^a = 3 : 3542 


iV e = 7, iVi x iV 2 = 7 x 7 


N A = 2 : 882 


N A = 2 : 882 


N A = 3 : 7105 


N A = 3 : 7105 


bosons atf = l/4inC = 3 band 


system size (Nb, Ni x N 2 ) 


•A/pis (-^4 , A^i , A^ 2 ) 


J^(N A ,N^ = N 1 N 2 ) 


N b = 4, Ni x % = 4 x 4 


A^a = 2 : 72 


N A = 2:72 


JVi, = 5, TVi x iV 2 = 4 x 5 


Aa = 2 : 130 


N A = 2: 130 


N b - 6, JVi x N 2 - 4 x 6 


A^a = 2 : 204 


N A = 2 : 204 


A^a = 3 : 680 


Aa = 3 : 728 


JVi, - 7, JVi x N 2 - 4 x 7 


Aa = 2 : 294 


Aa = 2 : 294 


N A = 3: 1428 


N A = 3: 1428 


A?6 - 8, Ni x iV 2 - 4 x 8 


N A = 2 : 400 


N A =2: 400 


N A = 3 : 2464 


N A = 3: 2464 


Af A = 4 : 7088 


A^a = 4 : 7752 


bosons at ^ = 1/5 in C = 4 band 


system size (Nb, Ni x N 2 ) 


Af£ s (N A ,N u N 2 ) 


M${N A ,N <t> = N 1 N 2 ) 


N b = 4, JVi x N 2 = 4 x 5 


N A = 2 : 90 


N A =2: 110 


iV 6 = 5, JVi x N 2 = 5 x 5 


A^a = 2 : 200 


A>a = 2 : 200 


Nb = 6, JVi x A?2 = 6 x 5 


Aa = 2 : 315 


Aa = 2 : 315 


A?a = 3 : 1000 


N A = 3 : 1360 


N b - 7, JVi x N 2 - 7 x 5 


N A = 2 : 455 


N A = 2 : 455 


AT A = 3 : 2345 


A?a = 3 : 2695 



